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ABSTRACT 

We present the results of X-ray spectral analysis of the superbubble around 
the OB association LH9 in the H II complex Nil in the Large Magellanic Cloud. 
Using observations from Suzaku, we have modeled the X-ray emission over the 
energy range 0.2-7.5 keV. We contrained the thermal spectrum below 2 keV using 
a thermal plasma model found in a previous XMM-Newton study of the Nil 
region. We find that the hard X-ray emission (> 2 keV) requires a nonthermal 
power-law component. The photon index of this component was found to be 
r = 1.72 ± 0.15. We performed an energy budget analysis for Nil using the 
known stellar content of LH9. We found that the observed thermal and kinetic 
energy in the superbubble is only half of the expected mechanical energy injected 
by stars. 

Subject headings: galaxies:ISM - ISM: Large Magellanic Cloud - X-rays: ISM 



1. Introduction 

Superbubbles (SBs) are formed by the combined actions of stellar wi nds from massive 
stars in OB associations and the eventual supernovae (SNe) of those stars (IBruhweiler et al. 
1980l ). The physical structure of a SB is very s imilar to that of a b ubbl e blown by the stella r 
wind of an isolated ma ssive star, as modele d by lCastor et al.l (Il975l ) and lWeaver et al.l (Il977l ). 
and more recently by iPittard et al.l (120011 ). The interior of a SB is shock- heated by stellar 
winds and SN ejecta to ~ 10^ — 10^ K and therefore exhibits diffuse X-ray emission. These 
regions of rarefied hot gas reach diameters of ~ 100 pc. 



^Current Address: Columbus State University, Coca-Cola Space Science Center, Columbus, GA 31907- 
5645 
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The study of SBs in the Galaxy is impeded by the effects of obscuration and confusion 
in the Galactic plane, making faint, extended objects such as SBs difficult to observe, partic- 
ularly at X-ray wavelengths. In contrast, the high Galactic latitude and low internal column 
density of the Large Magellanic Cloud (LMC) allows us a relatively un-obscured view of its 



popu lation of SBs, and as these SBs are all at roughly equal, known distances (50 kpc; [Feast 
19991 ). we can translate our observations into the physical properties of the SBs. 



Nil is the second largest H II region in t he LMC after 30 Doradus , and may constitute 
a more evolved version of this latter nebula (jWalborn fc Parked Il992l). It harb ors several 
associations of massive stars: LH9, LHIO, LH13 and LH14 (ILucke &: Hodgelll970l ). Its struc- 
ture is complex and reflects the interactions between the stars and their environment. The 
combined action of stellar winds and SN explosions from the central cluster LH9 has carved a 
hollow cavity in the surrounding ISM and created a SB shell ~ 120 pc in diameter. The hot 
shocked winds/eject a that fill this cavity emit X-rays and provide the pressure to drive the 
SB shell expansion (IMac Low et al.l 119981 ). Only SNe outside the SB can p roduce distinct 



SN re mnants (SNRs) such as NHL at the western edge of the Nil complex (jWilliams et al. 
1999h . 



Previous studies of the Nil region with ROSA T and XMM-Newton have presented th e 
soft X-ray spectral properties of the SB around LH9 (IMac Low et al.lll998l : lNaze et al.ll2004j ). 
Each study concluded that the soft X-ray component was dominated by a thermal emission 
characterized by a temperature of kT ~ 0.2 keV. Investigations into the high-energy (> 2 
keV) properties of Nil were limited in both of these studies: ROSAT was not sensitive to 
photon energies higher than ~ 2 keV; the XMM-Newton observations suffered from high 
particle background at these higher energies. The recently launched Suzaku has a moderate 
angular resolution, ~ 2', and greater sensitivity to higher energy X-ray photons, providing 
an excellent tool to study the high-energy diffuse emission in SBs. 

In this paper we present the analysis of the X-ray spectrum of the SB around LH9 in Nil 
using new Suzaku observations. In Section 2 we will provide a description of the observations 
and data processing procedures. Section 3 will describe the spectral analysis of the main SB 
in Nil based on the Suzaku observations. In Section 4 we will use the spectral results to 
perform an energy budget analysis to compare the observed thermal and kinetic energies of 
the SB to the energy input based on the stellar population of the OB association LH9. We 
will also discuss the possible mechanisms for the oberved non-thermal X-ray emission. 
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Observations and Data Reduction 



The X-ray observatory Suzaku {ASTR0-E2; iMitsuda et al.l 120071) includes five X-ray 
telescopes (XRTs), sensitive to soft X-rays up to ~ 10 keV (ISerlemitsos et al.l 120071 ). At 
the foci of four of the XRTs (XRT-I) are charge -coupled devices (CCD), known as X-ray 
Imaging Spectrometers (XIS; iKoyama et al.l 120071 ). Each CCD chip has an imaging area of 
1024 X 1024 pixels, with a corresponding pixels size of 24/im x 24/im. Three XIS CCDs (0,2 
and 3) are front-illuminated (FI), and one (XISl) is back-illuminated (BI). The field-of-view 
for the XIS CCDs are 18' x 18' with effective areas of 340 cm^ (FI) and 390 cm^ (BI) at 1.5 
keV. Suzaku is located in a lower orbit than XMM-Newton and Cha ndra, providing a low er 
and more stable particle background than other X-ray instruments (IMitsuda et al.l 120071 ). 



Nil was observed with Suzaku on 2006 November 7-8 (ID 501091010; PI Williams) for 
a total of 66.5 ks with a nominal pointing of R.A. 4^56'^50!90, Dec -66°24'21'.'6 (J2000). 
Data were acquired on the three front-illuminated imaging CCDs (XISO, XIS2 and XIS3) 
and the one back-illuminated CCD (XISl). The data were processed with the HEASoft 
v6.6 software suite. We first applied standard filters for bad pixels and event grade^. We 
removed time intervals when the source elevation was less than 10° above the Earth's limb. 
Additionally, we removed times for when the satellite was above or within 436s of the South 
Atlantic Anomaly. Finally, we screened out events from the onboard ^^Fe calibrator sources. 
After all of the data screening, we obtained a final exposure time of 22.2 ks. 

For the XIS instrumental response, we locally generated response matrices and ancillary 
response files using the tasks xisrmfgen and xissimarfgen, respectively. The XIS background 
spectra were taken from a blank sky observation toward the north ecliptic pole region, 
conducted for 95 ks on 2005 September 2-4. Non-X-ray background spectra were taken 
from night earth (NTE) observations compiled and available through the Suzaku Science 
Team. The NTE spectra were weighted according to the magnetic cut-off rigidity of the 
source observation. Source spectra were extracted using XSELECT from an elliptical region 
encircling the whole expanding SB around LH9. The spectra were then binned with the task 
grppha to include a minimum of 10 counts per bin over the entire energy range. 

An XMM-Newton observation of Nil was retrieved from the XMM-Newton data archive. 



Analysis of this observation was previously p r esente d in iNaze et al.l (120041 ). Our data pro- 
cessing procedure was identical to lNaze et al.l (120041 ). and our spectral analysis results were 
consistent with theirs. We used the resulting thermal model to constrain the low energy 
spectrum of the Suzaku observations. 



^As described in Suzaku ABC Guide: jhttp://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 
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Figure [T] shows the distribution of soft (0.2-l.OkeV) X-ray emission overlaid on an Ha 
image of the Nil complex from the Magellanic Cloud Emission Line Survey (MCELS; 



Smith et al.l Il998l ) . The XIS contours were generated using a combined image that was 



binned by a factor of 8 and smoothed with a Gaussian to the effective resolution of the 
Suzaku telescope ~ 2'. The peak soft X-ray emission in the Suzaku and XMM-Newton im- 
ages are offset, but lie within the ~ 30" pointing accuracy of Suzaku. The emission peak lies 
within the evacuated shell of LH9, where the stellar winds of the young stars have blown 
away the natal gas that formed the cluster. The three strong point sources seen in hard 
X-rays (Figure 2) lie outside our spectral extraction region. 



Results 



Due to the difference in angular resolution between Suzaku and XMM-Newton, 2' and 
~ 4" — 6", respectively, our extraction regions needed to be quite large. We extracted raw 
spectra from elliptical regions containing the entire SB around the OB association LH9. The 
region measured 9' x 7' with its major axis aligned along the east-west direction. 

The background subtraction in the XMM-Newt on spectra yielded count rates too low 
for meaningful fits above 2 keV, as was the case in iNaze et al.l ( 120041 ). For these spectra, 
we limited our fitting to the energy range 0.4-2 keV. The XIS spectra yielded a different 
challenge. Above 8 keV, the subtraction of the non- X-ray background was problematic, with 
incomplete subtraction of Ni Ka and Kf3 fluorescence lines. This was not alleviated with 
the subtraction of the exposure corrected blank-sky data. For these reasons, we limited our 
spectra fitting range to 0.2-7.5 keV. 

We fit the XMM-Newton and Suzaku spectra separately using XSPEC vl2.5.0 dis- 
tributed with HEASoft 6.6. The XMM-Newton data were used to contrain the low energy 
{E < 2 keV) thermal spectrum of Nil. We then used the resulting thermal model as a fixed 
input to the Suzaku model fi t. All fits include a fixed absorption column of 4.25 x 10^° cm~^ 
in the direction of the LMC (jPickey fc Lockmanlll990l ). An additional absorption component 
was used as a free parameter to account for the l ocal LMC column at the position of Nil. 
The first model fit attempted was a single mekal (IMewe et al.lll985r ) ther mal plasma model, 
the model found best for the XMM-Newton spectra in iNaze et al.l (120041 ). This single tem- 
perature model was unable to fit the Suzaku spectra above 2 keV. We then added a second 
thermal component using various models in order to fit the high-energy region of the spectra. 
An additional thermal component with a temperature of kT = 0.8 keV was able to improve 
the fit at low energies {E < 2 keV), but above 2 keV the models diverged quickly from the 
observed high-energy spectrum. The only model that produced an improved fit, the non- 
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equilibrium ionization model nei, yielded a gas temperature that was unphysical [kT > 15 
keV) for a system of this age. We therefore ruled out a second thermal contribution to the 
spectra. 

The best fit was achieved with a two-component model consisting of a thermal plasma 
(mekal) and non-thermal power-law. Initially, we allowed the metallicity to vary freely. The 
resulting fits did not deviate significantly from the average LMC abundances of 0.3 Zq. This 
is consistent with our understanding of this system as containing only a small fraction of SN 
enriched material intermix ed wit h the t ypical LMC ISM. This also agrees with the results of 
self-enrichment models by iNava I (120081 ) . who found that for OB associations with LMC-like 
initial conditions no measurable change in metal abundances can be seen in 4-5 Myr, the 
age of LH9. These current data are not sufficient to constrain abundance variations through 
SN enrichment over the broad spatial regions required by the Suzaku spatial resolution. We 
there fore fixed the abundanc e parameter to 0.3 Zq for the final model, typical of the LMC 
ISM faussell fc Dopitalll992h . The total absorbed fiux in the model is (5.5 ± 0.8) x lO^^^ 



erg s ^ cm ^. 

The model indicates that the thermal component of the gas has a temperature of kT = 
0.18 ±0.07 keV, which is in good ag reement with th e prop erties of other LMC superbubbles. 
The error bars are 1-a uncertainty. iMac Low et al.l (119981 ) obtained a temperature range of 
0. 1-0.2 keV to 99% confidence in the i?(95'^r observation. The temperature determination 
of iNaze et aP J2004h was kT 



0.181°:}^ keV. 



At higher energies (> 2 keV), a non-thermal component is needed to fit the X-ray 
spectrum. The photon-index {E oc i^^^) of the best-fit model is F = 1.72 ± 0.15, harder 
than would be expected for a synchrotron mechanism. The fiux of the nonthermal emission 
is (3.0 ± 0.8) X 10~^^ erg s~^ cm~^, about 35% of the total fiux. Figure [3] shows the final 
spectral model with the Suzaku and XMM-Newton data. The best-fit parameters are shown 
in Table 1. 



4. Discussion 

We have used the observed stellar content of the OB association LII9 in the Nil SB to 
estimate the mechanical energy injected by stars via fast stellar winds and SN explosions. 
We have also used optical and radio observations to determine the kinetic energy of the dense 
ionized and neutral gas in the SB, and X-ray observations to assess the thermal energy of 
the hot interior gas. The comparison of these quantities were then used to determine an 
energy budget for the SB. Finally we discuss the observed nonthermal X-ray emission and 
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possible physical mechanisms for this emission. 



4.1. Stellar Energy Input 



To estimate the stellar energy contri b utioii for Nil, we have considered the population 
of stars in LH9 as given by iParker et al.l (jl992l ). In that study, bolometric magnitudes and 
effective temperatures (as determined by their spectral types or UBV colors) were used 
to determine stellar masses and the present-day mass function (PDMF), then the PDMF is 



furthe r used to determine the initial mass function (IMF) of the OB association. iParker et al. 
(I1992I ) noted completeness problems for stars with M < 8 Mq. Ignoring the highest-mass 



bin due to poten tial evolutiona r y effe cts, iParker et al.l found the IMF to have a slope a = 
-1.45±0.1 while iMassev et all fll995h found a IMF slope of a = -1.38±0.16. We have 
adopted a middle value for the IMF of a = —1.4. 



We have used the Starburst99 v5.1 package (ILeitherer et al.lll999l ) to model the stellar 
energy inputs of Nil. Starburst99 requires three basic inputs: an IMF slope, the total mass 
of stars formed by the OB association in the 1-100 Mq range (essentially a scaling factor for 
the IMF), and a timescale to integrate the stellar energy inputs over. We use our adopted 
IMF slope of a = —1.4, and the observed main sequence star counts, to determine the total 
mass of stars formed by the OB association LH9. We note that the IMF will take the form 
/(M) = KM"""^ = ifM~^'^, where is a constant that can be determined from the star 
counts. Using the star counts in the 12-15Mq, 15-20, 20-25, and 25-40 Mq bins from 
Parker et al.l (119921 ). we find K =3500±200. The error bar reflects the uncertainty caused 
by integrating over different subsets of the total mass range. To find the total mass of stars 
formed by the OB association, we integrated f{M)dM over the 1-100 Mq mass range and 
find a total star mass of 7400±400 Mq. 

To estimate a timescale, we note that the earliest-type star observed in LH9 is an 
06 V star, with a mass of 40-6 Mq. Such a star would have a main sequence lifetime of 
~4-5 Myr (ISchaerer et al.lll993l ). and - assuming a single burst of star formation - gives us 
an estimate of the age of LH 9. The Starburst99 model predicts a total wind input energy 
of 8.6±0.5xl0^^ ergs at 5 Myr and that 8±1 supernovae have already occurred, adding 
8±lxl0^^ ergs of mechanical energy. The total mechanical energy from LH9, as predicted 
from the star counts, is therefore 17±2xl0^^ ergs. We assume that all of this energy has 
gone into forming the SB. 
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4.2. Observed Thermal and Kinetic Energy 

Using the best-fit model parameters, we are able to calculate some of the physical condi- 
tions within the SB. Since we used the XMM-Newton observations to constrain the thermal 
spectral fits, we have included those data in our calculations, leading to physical ranges. It 
should be noted that since our extraction region is very large, the physical estimates will 
contain hot gas that is not associated with LH9. The figures given here should be regarded 
as upper limits to the physical parameters within the main shell of Nil. 

We can estimate the electron density and mass of the thermal gas from the normahzation 
factor A for the spectral fits: 

where D is the distance to Nil, is the electron density, riu is the hydrogen density, and 
V is the volume in cgs units. We assume that the hydrogen and helium is fully ionized, and 
that riHc/^H = 0.1; giving us a total particle density of 1.92ne. The electron density and hot 
gas mass are then calculated by the expressions 

rie = 3.89 X 10^L»lV2y-i/2y-i/2 ^^-3^ 

Mgas = l.llnemnVf g, 

where / is the volume filling factor. From these quantities, along with the plasma temper- 
ature, we can derive the thermal energy and pressure of the hot gas around LH 9 by the 
following expressions: 

Etf, = 4.60 X 10-V(A;T)y/ erg, 
Pth = 3.05 X IQ-^eikT) dyne cm"^. 



We calculate that the X-ray emitting gas has an electron density of = 0.098 ± 
0.003 /"^'^^ cm~^, and that the total thermal energy for the gas is 2.3±0.9 x 10^^/^/^ ergs. The 
hot gas within the extraction region is calculated to have a total mass 

of 2762 ±85/1/2 

Finally, the thermal pressure of the emitting gas was calculated to be 5.38±2.14 x lO"^^/"^/^ 
dyne cm"^. 

It is most likely that the volume filling factor / will range in value from 0.5 — 1.0 due 
to the low density environment of the SB interior. If we adopt a value of / = 0.75 for the 
X-ray emitting gas, a medium value for the likely range, we calculate that the total observed 
thermal energy in the SB is 2.0 ± 0.8 x 10^^ erg. 
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Rosado et al.l (119961 ) found that the Nil SB had an Ha e xpansion velocity of 45 km 



in the shelL For a shell mass of 8.3 x 10^ Mq (jPunnd 120071 ). we calculate that the H II 
kineti c energy cont ribution is 1.7 x 10^^ erg. A survey of the H I distribution of SBs in the 
LMC (jDunndl2007l ) found that the neutral gas bubble was expanding at a rate of 25 km s~^. 
With an H I mass of 4.5 x 10^ Mq we calculate a kinetic energy of 2.8 x 10^^ erg for the 
neutral shell. So the total kinetic energy of the ionized and thermal gas is 4.5 x 10^^ erg, 
and the total observed energy (thermal plus kinetic) of the SB is 6.5 ± 2.5 x 10^^ erg. 



4.3. Energy Budget 



The calculated thermal and kinetic energy (6.5 ±2.5 x 10^^ erg) of the Nil SB represents 
an energy deficit of about 62 ± 13% when compared to the mechanical energy injected by 
the underlyi ng OB associa t ion L H9 (17 ± 2 x 10^^ erg). This deficit is similar to that found 
for N51D by lCooper et al.l (120041 ). In that case, it was shown that the observed thermal and 
kinetic energy were about 3 times lower than the stellar energy input. Suggested energy 
loss mechanisms were a large scale blow out, where the hot gas breaks through the shell 
and escapes into the ambient ISM, and evaporation of dense cool gas into the hot interior 
of the SB. Neither N51D nor Nil show evidence for blowouts, and radiative losses due to 
the interface of cool and hot gas at the shell are not sufficient to make up the energy deficit 
JCooper et al.ll2004h . 



4.4. Nonthermal X-ray Emission 



Our analysis of Suzaku observations of the SB Nil has demonstrated the presence of 
a nonthermal high-energy component to the spectrum that was not detected in previous 
ROSAT and XMM-Newton observations. This adds Nil to a small list of regions around 
OB associations th a t exhibit diffuse nont hermal X-ray emission: the SBs 30 Doradu s C 
jBamba et all booi . ISmith fc Wa"^ l2004h a nd DEM L192 (N5i p. ICooper et all boO^ in 
the LMC, and ROW 38 JWolk et all bom and Westerlund 1 JMuuo et all l2006h in the 
Milky Way galaxy. Proposed explanations for the nonthermal emission in these objects 
are synchrotron radiation or inverse Compton scattering of starlight. Each requires the 
acceleration of electrons to relativistic energies. Repeated SN shocks, turbulence in SNRs 
and colliding s tellar winds, all exp ected in OB associations, can accelerate the electrons to 
these energies (jParizot et al. 2004j ). 



In the case of Nil, the photon index (~ 1.72) is harder than would be expected for 
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a aynchrtron origin. It is more consistent with an inverse compton mechanism. The low 
angular resolution of the Suzaku observation makes it difficult to separate the individual 
contributions to the power-law emission, such as emission contributed by the stellar content 
of LH through the interaction of the stellar winds with the surrounding H II region, or SNR 
shock interactions with the SB shell. The raw XMM-Newton data indicate the presence of 
about six faint point sources within the spectral extraction region. Four of these sources 
seem to have relatively hard X-ray colors, but upon subtraction of background events, there 
are insufficient remaining photons to do spectral analysis on the sources. At least two of the 
XMM-Newton point sources contribute to the 1.0-2.0 keV emission along the NW rim of the 
SB. In the same region, one point source appears to contribute a few counts at energies >2 
keV in the Suzaku images. The point sources contribute less than ~ 1% of the total hard 
X-ray flux that we observe in the SB X-ray emission, an insignificant fraction of the total. 

The detection of nonthermal emission in Nil, combined with previous SBs, demon- 
strates that faint nonthermal emission may be present in many SBs where it has been 
undetected before due to limited high- energy sensitivity. The discovery of a high-energy 
nonthermal component in Nil indicates that such components may not be uncommon in 
SBs; further observations with instruments such as Suzaku may be expected to increase the 
sample of such objects. The origins of such components are not yet fully understood - and, 
indeed, multiple mechanisms may be at work. If in general they are shown to be due to the 
presence of recent internal SNR shocks, as is indicated for Nil, this will have interesting 
implications for the study of super bubble evolution and their role as cosmic-ray sources. Due 
to the complex nature of SB environments, it would be expected that nonthermal emission 
has multiple origins. With a larger sample we will be able to study more closely the question 
of what mechanism are generating this emission. 



5. Summary 

We have obtained Suzaku observations of the H II complex Nil in the LMC. Prominent 
diffuse X-ray emission is detected from the central SB. We have analyzed the Suzaku XIS 
data in conjunction with the XMM-Newton EPIC data. Wc find a main thermal emission 
component with a temperature of kT ~ 0.18 keV, consistent with previous analyses of 
RO SAT and XMM-Newton observations alone. Owing to Suzaku^s high sensitivity at photon 
energies > 3 keV, we arc able to extend the spectral analysis to hard X-rays and find a 
nonthermal X-ray component with a photon index of F ~ 1.72. This is consistent with 
an inverse compton mechanism, though we cannot rule out other mechanisms at this time. 
This nonthermal component contributes ~ 35% of the total X-ray emission from the analyzed 
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region. 

We have performed an energy budget analysis for Nil using the known stellar content 
of LH9. The total stellar wind energy and supernova explosion energy injected by LH 9 into 
the SB in the past 5 Myr are 8.6±0.5xl0^^ ergs and 8±lxl0^^ ergs, respectively. The total 
thermal and kinetic energies retained in the Nil SB are 2.3± 0.9 xlO^^ erg and 4.5x10^^ erg, 
respectively. We find a deficit of ~ 62% in the thermal and kinetic energy stored in the 
superbubble from the expected mechanical energy injected by the stars. This is similar to 
that observed in N51D. This could mean that these energy deficits may be common in SBs. 
More sensitive observations of SBs will help answer this question. 

The authors would like to thank the anonymous referee for helping to improve the 
quality of this paper. This work was supported in part by the following grants: Suzaku grant 
NNX07AF61G, LTSA grant NNG05GC97G (RMW). 

REFERENCES 

Bamba, A., Ueno, M., Nakajima, H., & Koyama, K. 2004, ApJ, 602, 257 
Bruhweiler, F. C., Gull, T. R., Kafatos, M., & Sofia, S. 1980, ApJ, 238, L27 
Castor, J., McCray, R., & Weaver, R. 1975, ApJ, 200, L107 
Chu, Y.-H., & Mac Low, M.-M. 1990, ApJ, 365, 510 

Cooper, R. L., Guerrero, M. A., Chu, Y.-H., Chen, C.-H. R., k Dunne, B. C. 2004, ApJ, 
605, 751 

Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 

Dunne, B. C. 2007, PhD thesis. University of Illinois at Urbana-Champaign 

Feast, M. 1999, in lAU Symposium, Vol. 190, New Views of the Magellanic Clouds, ed. Y.-H. 
Chu, N. Suntzeff, J. Hesser, & D. Bohlender, 542 

Koyama, K., et al. 2007, PAS J, 59, 23 

Leitherer, C, Schaerer, D., Goldader, J. D., Delgado, R. M. G., Robert, C, Kune, D. F., de 
Mello, D. F., Devost, D., & Heckman, T. M. 1999, ApJS, 123, 3 

Lucke, P. B., & Hodge, P. W. 1970, AJ, 75, 171 



- 11 - 



Mac Low, M.-M., Chang, T. H., Chu, Y.-H., Points, S. D., Smith, R. C, & Wakker, B. P. 

1998, ApJ, 493, 260 

Massey, P., Lang, C. C, Degioia-Eastwood, K., & Garmany, C. D. 1995, ApJ, 438, 188 
Mewe, R., Gronenschild, E. H. B. M., & van den Oord, G. H. J. 1985, A&AS, 62, 197 
Mitsuda, K., et al. 2007, PAS J, 59, 1 

Muno, M. P., Law, C, Clark, J. S., Dougherty, S. M., de Grijs, R., Portegies Zwart, S., & 
Yusef-Zadeh, F. 2006, ApJ, 650, 203 

Nava, A. H. 2008, Ph.D. Thesis, U of Oklahoma 

Naze, Y., Antokhin, I. I., Rauw, G., Chu, Y.-H., Cosset, E., & Vreux, J.-M. 2004, A&A, 
418, 841 

Parizot, E., Marcowith, A., van der Swaluw, E., Bykov, A. M., & Tatischeff, V. 2004, A&A, 
424, 747 

Parker, J. W., Garmany, C. D., Massey, P., & Walborn, N. R. 1992, AJ, 103, 1205 

Pittard, J. M., Dyson, J. E., & Hartquist, T. W. 2001, A&A, 367, 1000 

Rosado, M., Laval. A., Le Coarer, E., Ceorgelin, Y. P., Amram, P., Marcelin, M., Goldes, 
G., & Gach, J. L. 1996, A&A, 308, 588 

Russell, S. C, & Dopita, M. A. 1992, ApJ, 384, 508 

Schaerer, D., Meynet, G., Maeder, A., & Schaller, G. 1993, A&AS, 98, 523 
Serlemitsos, P. J., et al. 2007, PAS J, 59, 9 
Smith, D. A., & Wang, Q. D. 2004, ApJ, 611, 881 

Smith, R. C, & MCELS Team 1999, in lAU Symposium, Vol. 190, New Views of the 
Magellanic Clouds, ed. Y.-H. Chu, N. Suntzeff, J. Hesser, & D. Bohlender, 28 

Walborn, N. R., & Parker, J. W. 1992, ApJ, 399, L87 

Weaver, R., McCray, R., Castor, J., Shapiro, P., & Moore, R. 1977, ApJ, 218, 377 

WiUiams, R. M., Chu, Y.-H., Dickel, J. R., Smith, R. C, Milne, D. K., & Winkler, P. F. 

1999, ApJ, 514, 798 



- 12 - 

Wolk, S. J., Spitzbart, B. D., Bourke, T. L., & Alves, J. 2006, AJ, 132, 1100 



This preprint was prepared with the A AS IM^jX macros v5.2. 



Table 1. 



Model Parameter LH9 



Column Density Nh (10^^ cm ^) 






0.52 ±0.04 


Plasma Temperature kT (keV) 






0.18 ±0.07 


Thermal Flux (absorbed) (lO"^^ 


erg s 


cm ^) 


25.4 


Thermal Flux (unabsorbed) (10" 


-14 

erg 


; s~-^ cm~^) 


954 


Photon Index (7) 






1.72 ±0.15 


Nonthermal Flux (absorbed) (10 


— 14 

erj 


g s~^ cm~^) 


29.6 


Nonthermal Flux (unabsorbed) ( 




erg s~^ cm~^) 


46.8 


(d.o.f) 






544 (442) 
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Fig. 1. — [left) Suzaku soft (0.2-1.0 keV) X-ray contours overlaid on an Ha image of Nil. 
The locations of the OB associations LH9, LHIO and LH13 are indicated by the arrows. The 
peak of the X-ray emission is coincident with the position of LH9. {right) XMM-Newton soft 
X-ray contours for the same region. The SNR NllL to the right of the SB lies outside the 
Suzaku field-of-view. 
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Fig. 2. — Combined XIS images of Nil with X-ray contours overlayed. Energy ranges for 
each image are (a) 0.2-l.OkeV, (b) 1.0-2.0keV and (c) 2.0-lO.OkeV. 
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Fig. 3. — (top) Combined spectra from all XMM-Newton and Suzaku detectors. Colors 
denote the spectra from each instrument. Solid lines show best-fit model, {bottom) Plot of 
the Suzaku spectra with the contributions from each model. The spectral features below 2 
keV represent emission line blends of highly ionized iron, neon and magnesium. 



